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We study analytically the dynamic behaviors of quantum correlation measured by quantum discord between 
two uncoupled qubits, which are immersed in a common Ohmic environment. We show that the quantum discord 
of the two noninteracting qubits can be greatly amplified or protected for certain initially prepared X-type states 
in the time evolution. Especially, it is found that there does exist the stable amplification of the quantum discord 
for the case of two identical qubits, and the quantum discord can be protected for the case of two different qubits 
with a large detuning. It is also indicated that in general there does exist a sudden change of the quantum discord 
in the time evolution at a critic time point t^, and the discord amplification and protection may occur only in the 
time interval < ? < f^ for certain X-type states. This sheds new light on the creation and protection of quantum 
correlation. 
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I. INTRODUCTION 

It is well known that the total correlation in a bipartite quan- 
tum system can be measured by quantum mutual informa- 
tion fr, Vl, which may be divided into classical and quan- 
tum parts [2-9]. The quantum part is called quantum discord 
which is originally introduced by Olliver and Zurek |l7|]. Re- 
cently, it has been aware of the fact that quantum discord is 
a more general concept to measure quantum correlation than 
quantum entanglement since there is a nonzero quantum dis- 
cord in some separable mixed states [7J. In fact, quantum 
discord is a different type of quantum correlation than en- 
tanglement, and it can be considered as a more universal re- 
source than quantum entanglement in some sense. As shown 
in Refs. lliollllll . although there is no quantum entanglement, 
quantum discord can also be responsible for the quantum com- 
putational efficiency of deterministic quantum computation 
with one pure qubit lll2ll . In addition, much recent attention 
has been paid to many relative topics of quantum discord Iil3r - 
I22II . such as quantum discord of open quantum systems. 

We know that any realistic quantum systems interact in- 
evitably with their suiTounding environments, which intro- 
duce quantum noise into the systems. As a result, the quan- 
tum systems will lose their energy (dissipation) and/or coher- 
ence (dephasing). Thus it is of fundamental importance to 
know the influence of the environment on quantum correla- 
tion. In several recent papers ll23l - l27ll . the quantum correla- 
tion dynamics in open quantum systems have been studied. It 
was shown that the quantum correlation measured by quantum 
discord is more resistant against the environment than quan- 
tum entanglement |23]. For a certain class of states under 
Markovian dynamics, the quantum entanglement can disap- 
pear within a finite time, a phenomenon called by entangle- 
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ment sudden death f28'l, which has been widely investigated 
in recent years [28. - 3 3.1 . Differently, quantum discord only 
vanishes asymptotically at infinite time ll23ll34|[ . Moreover, 
for some special initial states, quantum correlation in a bipar- 
tite quantum system will not be affected by the decoherence 
environment during an initial time interval jTM . 

Based on the above mentioned two facts that quantum dis- 
cord is a useful resource for quantum information processing, 
and that quantum systems couple inevitably with their envi- 
ronment. Naturally there are two interesting questions: (i) 
if the environment can enhance the quantum discord of the 
systems to realize the discord amplification? (ii) if we can 
obtain a stable quantum discord induced by the environment? 
With these questions, in this paper, we study the dynamics 
of quantum discord between two non-interacting qubits im- 
mersed in a common Ohmic environment. We show that it 
is possible to amplify and protect the quantum discord by the 
qubit-environment interaction under some conditions. Espe- 
cially, when the two qubits are identical, we find that the phase 
decoherence can induce a stable amplification of the initially- 
prepared quantum discord for certain X-type states. 

This paper is organized as follows. In Sec. [Ill we present 
our physical model and its solution. In Sec. Hill we investigate 
dynamical behaviors of quantum discord for the so-called X- 
type initial states. In Sec. IIVI we study analytically and nu- 
merically the quantum discord for certain X-type initial states, 
and show the initially prepared discord can be amplified or 
protected for X-type states. Finally, we conclude this work in 
Sec.|V] 



II. PHYSICAL MODEL AND SOLUTION 

Let us start with introducing the physical system, as shown 
in Fig.[T] two uncoupled qubits, qubit A and qubit B with en- 
ergy separations loa and ojb respectively, are immersed in a 
common environment. The Hamiltonian of the total system 
including the two qubits and the environment is composed of 




FIG. 1: (Color online) Schematic of our physical system: two un- 
coupled qubits, of energy separations oi^i and aig, are immersed in a 
common heat bath with temperature T . 



four parts, 



Hj — Hs + ti£ + Hi + Hf> 



(1) 



where Hs is the Hamiltonian of the two qubits. He is the 
Hamiltonian of the environment, Hj represents the interaction 
Hamiltonian between the two qubits and the environment, and 
Hr is a renormalization term, which is introduced originally 
in Ref. 1.36.1. 

The expression of the Hamiltonian Hs reads 



1 _ 1 



B' 



(2) 



where (T-^^^j 



\0)a(B) <0| - 1 1 )a(B) i 1 1 with |0>^(B) and 1 1 }^^s) be- 
ing the excited and ground states of the qubit A (B). Hereafter, 
we take ft - I. 

The environment of the two qubits is modelled by a heat 
bath with temperature T, which is composed of an infinite set 
of harmonic oscillators with the following Hamiltonian, 



He ^Yj ^kblbk. 



(3) 



where Uk is the frequency of the k\h harmonic oscillator de- 
picted by the usual bosonic creation and annihilation operators 
b'l^ and bk, satisfying the commutative relation {bk, /j^] = 1. 

As for the interaction Hamiltonian /// between the two 
qubits and the heat bath, we assume it is given by the fol- 
lowing expression 



Hi^HsYjgkibl+bk\ 



(4) 



where gk is the coupling constant between the system and the 
kxh harmonic oscillator of the heat bath. It is obvious that Hi 
commutes with the Hamiltonian Hs ■ Therefore there is no en- 
ergy dissipation in this system. Note that this form of coupling 
has been used to study quantum decoherence in Bose-Einstein 
condensation and trapped ion by Kuang and coworkers ||37I1 . 
The renormalization term has the following form 



H^^iHsfY, 



,2' 



(5) 



According to Ref. f37], Hamiltonian ([U can be exactly 
solved by making use of the unitary transformation 



U - exp 



HsJ]8k(bl-h) 



(6) 



Corresponding to Hamiltonian (HJ, the total density opera- 
tor of the system plus the heat bath can be expressed as 



pj(t) = e-'""U-'e 



iHs I fj- 1 -it Ei: <^kblbk 



''"• f/pr(0)f/"'e''2' t^thlbt f)g 



iH,t 



(7) 
where p7-(0) is the initial total density operator 

We assume the system and the heat bath are initially uncor- 
rected with the density operator p7-(0) =p(0) ® pr{Q), where 
p(0) is the initial density operator of the system, and Pr(0) 
is the density operator of the heat bath, which is assumed to 
be pr = Y\k ®pk{Q), where pk(Q) is the density operator of 
the A:th harmonic oscillator in thermal equilibrium of temper- 
ature T . We can obtain the reduced density operator of the 
system, denoted by pit) - TrR[p7-(f)]. Its matrix elements 
in the eigen representation of H^ (with the four eigenstates 
|00> , |01) , |10) , |1 1)) can be written as 

■^g-i{E{l',^J'„)-E{l/,,lBW ^ (g) 

where E{Ia, h) (£(1'^, I'g)) is the eigenvalue of the operator Hs 
with the corresponding eigenstate \Ia, Ib) d'^, I'g)), the expres- 
sion of EUaJb) is EUaJb) = [i-l)'''(^A + (-1)''>cjb]/2. The 
quantity Rg' / ){if,,ig){t) is a reservoir-dependent part given by 



R(V„)(lAjB)(t) 



-i[E'-(i'j')-EHiAJBmAi) 



xe 



-{E{l'J'B)-E{lA,lB)fQ2if) 



(9) 



The two reservoir-dependent functions Q\{t) and 22(0 in the 
above Eq. ^ are given by 



Qi{t) 



d(jjj{(jj) — T— sin(wO> 



A) W' 

22(0 = 2 f dojJioj)?-^ 
Jo w^ 



(10a) 



sin^(f)coth(f). (10b) 



Here we have taken the continuum limit of the reservoir 
modes Yjk ^ L dioJ{(jj), where 7(w) is the spectral den- 
sity of the reservoir, gico) is the corresponding continuum ex- 
pression for gk, and /3 - l/T (with the Boltzmann constant 

kB = 1). 



III. DYNAMICS OF QUANTUM DISCORD FOR 
TWO-QUBIT X-TYPE STATES 

In this section, we investigate dynamics of quantum discord 
for the two non-interacting qubits in the phase decoherence 
environment for the three-parameter two-qubit X-type states. 
Quantum discord [7] is defined as the difference between the 



total correlation and the classical correlation with the follow- 
ing expression 

D(p)^I(Pa:pb)-C{p). (11) 

Here the total correlation in a bipartite quantum state p is mea- 
sured by quantum mutual information given by 



I(Pa-Pb) = S(pa) + S (pb) -S (p). 



(12) 



where S (p) - -Tr(plogp) is the von Neumann entropy, 
pA - TrB(p) and ps - Tr/i(p) are the reduced density oper- 
ators for subsystems A and B, respectively. And the classical 
correlation between the two subsystems A and B can be de- 
fined as 



C(P) 



max 

{Pk] 



S(PA)-Y,PkS0A"> 



^ikh 



S (Pa) -min 

{Pt) 



Z 

L k 



PkSipf) 



(13) 



Here {Pk] is a set of projects performed locally on the sub- 
system B, and plf - ^Trg [(/^ ® Pk)p(lA ® Pk)] is the state 
of the subsystem A conditioned on the measurement of the 
outcome labelled by k , where p^ = TtabK^a ® Pk)p(jA ® Pk)} 
denotes the probability relating to the outcome k, and Ia de- 
notes the identity operator for the subsystem A. 

In terms of the relation given in Eq. ([8]), we can study the 
quantum discord dynamic properties of the two qubits. We 
assume the two qubits are initially prepared in a class of state 
with maximally mixed marginals {pa{B) - Ia{B)I'^) described 
by the three-parameter X-type density matrix 



P(0) 



( 3 

hB + 2_fi^A ® ^ 

i=l 

( \+Ci Cy-C2\ 

1 - C3 Ci H- C2 

Ci H- C2 1 — C3 

C\ - C2 1 -H C3 



(14) 



where Iab is the identity operator in the Hilbert space of the 
two qubits, &'^ and &'g (/ - 1,2,3 mean x,y,z correspond- 
ingly) are the Pauli operators of qubit A and qubit B, and c, 
(0 < \cj\ < 1) are real numbers satisfying the unit trace and 
positivity conditions of the density operator p. The density 
operator p includes the Werner states and the Bell states as 
two special cases. Under the decoherence environment, the 
evolution of density operator p(f) initially prepared in Eq. (fl4l i 
can be obtained according to Eq. dS). Its explicit form at time 
t is 



( \+Ci 



Pit) = 4 



\IJ.e 



-/A, 





1-C3 









1-C3 








1 -1-C3 



(15) 



where we have introduced the following parameters 



Ai 


= (WA + (^B)t, 


A2 = (oJa - OJB)t, 




/^(O 


= (Ci -C2)ri(0, 


V{t) = (Ci + C2)r2(0, 




ri(0 


_ ^-(tJA+bJBfQ2(t) 


y2(f) = e-(<^A~<^BfQ2it) _ 


(16) 



As shown in Eq. ( fTSI l. the decay parameters yi and 72 de- 
pend on the frequencies coa and cob of the two qubits and the 
reservoir-dependent function Q2. For convenience, we define 
a detuning parameter of the two qubits as 



Wa 
ODb' 



(17) 



which indicates that the detuning parameter r - \ for two 
identical qubits due to coa - cjb, and r ^ 1 for two different 
qubits due to oja + (jJb- In other words, the two qubits are 
resonant (non-resonant) when the detuning parameter r - \ 
(r ^ 1). Then the two decay parameters y\ and 72 can be con- 
nected with each other through the detuning parameter with 
the following simple expression 



72 = r 



iT^f 



(18) 



When the two qubits are identical, i.e., r - \, from Eq. (fTSI l 
we have A2 = 0, 72 - 1, and v - {c\ + C2). In this case, 
equation ( fTsl l indicates that there does exist a decoherence- 
free subspace with two basis states 1 1 , 0)ab and |0, 1 )ab for the 
X-type initial states under present consideration. 

In order to obtain the mutual information of state p{t) given 
in Eq. ( fTSl l, we first calculate the four eigenvalues of p(f), 

M,2 = -(I+C3+JU), ^3,4 = -(1 -C3+V). (19) 



Then the mutual information reads 



I(Pa-Pb) = 2 -h ^Ai log Ai. 



(20) 



(=1 



Note that here we have used S (pAit)) - S (pBit)) - 1 > since the 
two reduced density matrixes pA(t) and pB(t) are maximally 
mixed, that is Pa(0 - ^a/2 andpB(f) - Ib/2. 

For calculation of the amount for the classical correlation 
C(p) defined in Eq. ( fT3]) . we propose the complete set of or- 
thogonal projectors {Pk - \0k){0k\,k = ||, ±) for a local mea- 
surement performed on the subsystem B, where the two pro- 
jectors are defined in terms of the following two orthogonal 
states 



leii) = cos 6* |0) -He''* sin 011), 
|6»^) = e-'V-sin6»|0)-cos6»|l), 



(21a) 
(21b) 



with < 6* < k/2 and < (p < In. After the two project mea- 
surements with p„ - Pi_ - 112 , the reduced density matrices 
of subsystem A read 



Pa = 4 



1 /2(1-C3cos(20)) esin(26l) 

e*sin(26i) 2(1 -i- C3 cos(26l)) 



^i _ W 2(1 -H C3 cos(20)) -esin(26') 
^'^ ~ 4 1 -e* sin(20) 2(1 - C3 cos(20)) 



(22a) 
(22b) 

where we have introduced the parameter e = yue'**' '^^ -1- 
yg'(A2+0) According to Eq. ( l22l l. it is straightforward to obtain 
the eigenvalues of the reduced density matrix p^ as follows: 



^fl- 5(1 ± A), 



(23) 



where we have defined A as 



A 



cX cos 



'■(20) + -[i? +v^ + 2//vcos(A2 - Ai + 20)] 



xsin^iie)] 



(24) 



Since ^^ ^(\-^^, then we have ^CpJ^) = S(p^) = /(A) defined 
by 



/(A) = -^log/^-^ 



2 °^\ 2 / 2 
which leads to the classical correlation 

C(p(0) = l-min[/(A)]. 



1+A /1+A 

log2 



(25) 



(26) 



Since the function /(A) is a monotonically decreasing func- 
tion, therefore for getting the minimal value of S (p\) - /(A), 
we should choose proper parameters and </> to ensure the pa- 
rameter A defined in Eq. jTM is maximal. Obviously, from 
Eq. JTM we can see that the maximal value depends on C3,ju, 
and V. From Eq. ( l24l i it is easy to get the following inequality 



A < 



,2,„,2(2e)+M±M)^sin2(20) 



1/2 



|C3|, for |C3| > (\n\ + |v|)/2, 

(\H\ + |v|)/2, for |C3| < (\n\ + |v|)/2. 



If we define ;t'(f) as 

;^(0 = max[|c3|,(IM0l + |v(0l)/2] 
then the classical correlation can be expressed as 



CiPit)) = 2 



1 + i-iTx 



log2[l+(-l)";^]. 



(27) 



(28) 



(29) 



11=1 



Therefore, the quantum discord can be written as 



D(p(t)) = 2 + J^Ai\og2Ai-C(p(t)), 



(30) 



(=1 



where the amount of the classical correlation C{p(t)) is given 
by Eq. ( |29] |. In principle, we have obtained the dynamics of 
the quantum discord according the above expression given in 
Eq. ( |30] |, provided that we know the initial condition of the 
system. In what follows we will study dynamic properties of 
the quantum discord for some initial states in detail. 



IV. QUANTUM DISCORD AMPLIFICATION 

In the section, we would like to find the possibility of quan- 
tum discord enhancement induced by phase decoherence in 
the dynamic evolution of the X-type quantum states. From 
Eq. ( fT6l ) we can see that both /v(f) and v(f) are two decaying 
functions with respect to the evolution time t due to the reser- 
voir function 22(0 being a increase function indicated below. 



hence there may exist a critic time f^. At this time, the follow- 
ing equation may be satisfied 



Wc)\ + \v{t. 



= k3l, 



(31) 



which is the equation to determine t^. Obviously, the critic 
time tc depends on initial-state and reservoir parameters. 
From Eqs. ( fT6b and (l28b we can obtain 



xit) 



i \C3l 

\ i\ij(t)\ + |v(0l)/2. 



for t > tc, 
for t < tc. 



(32) 



which indicates that the classical correlation expressed by 
Eq. ( |29l ) exhibits a sudden change at the critic time tc deter- 
mined by Eq. ( [3TT i. The classical correlation decays mono- 
tonically from the initial time to the critic time tc, and then 
keeps constant after the critic time tc- Then making use of 
Eqs. (29-32) we can obtain expressions of the quantum dis- 
cord in different regimes of the time evolution 



4 2 

D (pit)) =2 + ^^,. log2 ^,- - ^ £ ^« ^°82 a;, t<t, 



2^ 

«=i 



(33a) 



1 
D (pit)) =2 + J^Ai log2 'i/ - 2 X ^" ^°S2 ^"' f ^ fc, 

i=l n=l 

(33b) 
where A„ and A'„ are given by 

A„ = 1 + (-1)" |C3| , a; = 1 + (-1)«M±M. (34) 



From Eqs. (I33ab and ( I33bb we can see that there does exist 
a critic time at which the quantum discord evolution exhibits 
a sudden change. Therefore, dynamic evolution of the quan- 
tum discord may have different behaviors in the time interval 
< f < fc and after the critic time. In what follows, we shall 
investigate dynamic characteristics of the quantum discord in 
the time evolution for two identical and diff'erent qubits, re- 
spectively. 



A. The case of two Identical qubits 

For the sake of simplicity, we first consider the case of two 
identical qubits and assume that the three initial-state param- 
eters obey C2 - and ci > cj, > 0. In this case, we have 
wa = wb = il, r = 1 and 72 = 1 ■ From Eq. (ISTT i we find that 
the critic time tc satisfies the equation 



^ 4Q2 \ c, 



(35) 



where Qiit) is the reservoir function defined in Eq. ( fTOb which 
is determined by the spectral density of the environment 7(w). 
We consider the Ohmic reservoir 13811 with the spectral density 
J(cl>) given by 



J(aj) 



rjo) 



expl 



(36) 
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FIG. 2: (Color online) Plot of the scaled critic time oi^fr given in 
Eq. (38) vs the qubit parameter r^Qr and the initial-state parameter 
C3/C1 at zero temperature. 



FIG. 3: (Color online) Plot of the quantum discord 2)„/(0 for the two 
identical qubits vs the scaled time 0)^1 and the initial-state parameter 
C[ at zero temperature when rjQr = 1. 



where w^ is the high-frequency cutoff, and 77 is a positive char- 
acteristic parameter of the reservoir. With this choice, at low 
temperature the function 22(0 is given by the following ex- 
pression 



Q2(t) = nl^\n[l+iojjf] 



In 



^sinhf^ 



TTt 



The the critic time tc can be obtained from Eqs. dSST l 
and (IJTT i. However, it should be pointed out that there does 
exist the critic time only for certain values of the three initial- 
state parameters (ci,C2,ct,), not for arbitrary values of these 
three parameters. Take into account the positivity of the den- 
sity operator and the reservoir function given by Eq. d37] i, we 
find that there exists the critic time te for some initial states 
with < ci/2 < C3 < ci < 2/3 and C2 = 0. It is a diflicult task 
to get an analytical expression of the critic time f^ for an ar- 
bitrary reservoir temperature. However, at zero temperature, 
from Eqs. (l35T l and (l37l i it is straightforward to get the analytic 
result of the critic time given by 



-V(2c3/c,-ir'/(2.n^-)-l. 



(38) 



From Eq. ( l38l l we can see that the critic time t^ depends on 
not only the characteristics of the two qubits and their environ- 
ment but also the initial-state parameters. For a set of given 
parameters (ci, C3, 77, Q), the lower is the cut-off frequency of 
the reservoir, the longer is the critic time. In particular, when 
2c3 - ci we have tc - 00. This means that the quantum dis- 
cord wiU not exhibit the sudden change in the whole process 
of the time evolution. In Fig.|2]we display the critic time with 
respect to the parameters rjQ? and C3/C1 at zero temperature. 
In the following we show that at zero temperature it is possi- 
ble to realize the discord amplification in the case of 2c3 - ci 
for the two identical qubits. 

For the two identical qubits, when 2c3 - ci and C2 = the 
critic time approaches the infinity. Then the quantum discord 



is given by 



4 ^ 2 

Didit) = 2 + ^ ^; log, ^i-:^J] K log2 a;,. < f < 



(39) 



(37) where Ai, A',, and ji if) are given by 



M.2 = ^[2 + ci(l+2ri(0)], ^3.4 = l[2-ci(l+2)], 



A'1.2 = i + 2'^i*^^ +ri(0), ri(0 



,-40^62(0 



(40) 



where the temperature-dependent reservoir function Q2{t) is 
given by Eq. ( [37l ). 

In Fig. [3] we have plotted the quantum discord Didif) with 
respect to the initial-state parameter c\ in the time evolution. 
Figure [3] clearly indicates the amplification of the initial dis- 
cord in the time evolution. And the discord amplification be- 
comes more apparent for large values of c\. In particular, 
when the time approaches to the infinity, we have ji - 0, 
hence the quantum discord reaches its maximal value given 
by 



£>,rf('») = ^-^log2(2 + c,)-^--^log2(2-ci) 



2-3ci 



log2(2 - 3ci), 



(41) 



which implies that the interaction between the reservoir and 
the two identical qubits can lead to the stable amplification 
of the initially prepared discord. And the stable value of the 
quantum discord depends on only the initial-state parameters. 
In Fig. ID we display the stable discord D„/(oo) with respect to 
the initial-state parameter c\. From Fig.|4]we can see that the 
larger is the initial-state parameter ci, the more apparent the 
amplification of the quantum discord. 

In order to calculate the amplification rate of the quantum 
discord, we need the quantum discord of the corresponding 





FIG. 4: (Color online) Plot of the quantum discord given in Eq. (41) 
with respect to the initial-state parameter C[ . 



FIG. 5: (Color online) Plot of the amplification rate of the quantum 
discord F = 2),j(oo)/D(0) with respect to the initial-state parameter 



initial states with the following expression 

D(0) = -l + ^^log2(2-ci)+^^log2(2 + ci) 

2 + 3ci 2 — 3ci 

+ —^ log2(2 + 3ci) + —^ log2(2 - 3ci) 

-i^log2(l+ci)-i^log2(l-ci). (42) 

Then we define the amplification rate of the quantum dis- 
cord as r = £),v/(oo)/D(0). In Fig. 5 we have plotted the 
amplification rate T with respect to the initial-state parameter 
ci . Figure|5]indicates that the initially-prepared discord under 
present consideration can be amplified for the whole regime 
in which ci can takes its values, and the amplification rate is 
larger in the middle regime than in other regimes for the values 
of ci. In particular, we can obtain the maximal amplification 
rate Fmax ~ 2.17 when the initial-state parameter ci ^ 3.3. 

Physically, the stable amplification phenomenon of the 
quantum discord for the two resonant qubits in the common 
heat bath is related to the existence of a decoherence-free sub- 
space of the two qubits in the time evolution. In fact, from 
Eq. (15) we can see that when the two qubits are resonant, 
i.e., r - 1, we have v(t) = 0, the subspace formed by |0, 1)ab 
and |0, 1 }ab is a decoherence-free subspace of the two qubits 
in the time evolution. Therefore, the stable amplification phe- 
nomenon of the quantum discord can be understood as a quan- 
tum coherent effect of the two qubits in the common heat bath. 



B. The case of two different qubits 

Dynamic evolution of the quantum discord for two uncou- 
pled different qubits in a common environment can be stud- 
ied through analyzing the first derivative of the quantum dis- 
cord with respect to time. In this case, we have r ?i 1 due to 
oja + (Ob- We now discuss the general conditions to realize 
the amplification of the quantum discord by phase decoher- 
ence. Note that the time dependence of the quantum discord 



given by Eq. ( |33] ) can be described in terms of only one pa- 
rameter y\ it) since the other time-dependent parameters /I,, 
A',, and j2{t) can be expressed by j\{t) through the relations 
given by Eqs. ( fT6l ) and ( fTTj l. For the Ohmic -reservoir case un- 
der our consideration and at low temperature, making use of 
Eq. ( l36b we can find ji (t) to be 



rt(0H[i-(-^0^l''£-hg 



-nix+xyiJ- 



(43) 



From Eq. (|44] |. it is easy to see that the first derivative of y\ 
with respect to time t is always negative, i.e.. 



dy\ 
dt 



<0, 



(44) 



which means that the parameter ji monotonically decreases 
with respect to time /. In fact, y\ - j2 - \ sX the beginning 
and ji — > 0, 72 — > at infinite time. 

Hence, the derivative of the quantum discord with respect 
to evolution time can be expressed as 



dD(p{t)) dD(p(t))d7i 



dt 



dyi dt 



which implies that when 

dp {pit)) 



<0, 



and making use of Eq. (45) we can find 
dD(p{t)) 



dt 



>0, 



(45) 



(46) 



(47) 



which indicates that the quantum discord is enhanced with the 
time evolution. This is the condition of the amplification of 
the quantum discord. In the following we will further study 
the amplification condition. 
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FIG. 6: (Color online) Plot of the quantum discord given in Eq. i33t 
vs the scaled time coj at zero temperature for different values of the 
detuning parameter r: r = I (solid black line), r = 1 .5 (dashed purple 
line), r = 2 (dotted blue line), and r = 5 (dot-dashed red line). Other 
parameters are set as ci = 0.6, C2 = 0, and c^ = 0.3. 



In order to demonstrate the possibility of the discord ampli- 
fication, we calculate the first derivative of quantum discord 
with respect to yi. Substitution of Eqs. ([18) , (O, and d34] i 
into Eq. ( |33] |. the first derivative of !D(p{yi)) with respect to 
yi has the following forms 



dy\ 



F(yi) + G(yi), 
G(ri), 



< t <tc, 

t>tc. 



(48) 



where we have introduced the following three functions 



F(7i) 
G(ri) 
Kyi) 



ci - C2\ + \h\ , i2-\n\' 

I0g2 



4 

C\ -C2 



V 



2 + |/i| + |v| 



log2 



1 + C3 + ;U 

1 -I- C3 - /I 



h 

4log2 



1 - C3 H- V 
1 - C3 - V 



r-Hl 



(49) 



For the initial state of < q < ci < 1 and |c3| < 1, from 
Eq. ( |49l ) we can see that the function G{y\) is always posi- 
tive. This implies that when t > tc the derivative of the quan- 
tum discord with respect to the time is always negative, i.e., 
dDipiyi)) Idt < 0. Hence, the quantum discord always de- 
creases monotonically when t > tc- That is, it is impossible 
to amplify the prepared quantum discord in the time evolution 
after the critic time tc- 

On the other hand, from Eqs. ( l48T l and ( |49l ) we can see that 
the function F{y\) is always negative due to (2 - \n\ - |v|)/(2 + 
M + |v|)< l.Then,5£)(p(ri))/5ri < 0, i.e., 5£)(p(ri))/5f < 
0, is possible due to the competing change of the two functions 
F{y\) and Giji) in the time evolution when t < tc- Therefore, 
it is possible to amplify the initially prepared quantum discord 
before the critic time. 

We now investigate numerically the discord amplification 
for certain initially prepared X-type states of two different 



qubits. In Fig. |6] we display the dynamic evolution of the 
quantum discord at zero temperature. Here the two qubits are 
initially prepared in the X-type state with c\ - 0.6, ci - 0, 
and C3 - 0.3, the detuning parameter r = 1, 1.5,2, and 5, 
respectively. The solid black line corresponds to the case of 
two identical qubits discussed in the previous subsection. In 
this case, the initially prepared discord can be amplified and 
this discord amplification is stable with the critic time tc - 00. 
From Fig. 6 we can see that the detuning of the two qubits 
described by the parameter r can affect seriously the discord 
evolution. The critic time at which the sudden change of the 
discord happens becomes shorten with the increase of the de- 
tuning parameter r. After the critic time, the discord begins 
to decay asymptotically. There exists an amplification regime 
of the discord before the critic time. The amplification regime 
becomes narrow with the increase of r, and the amplification 
disappears when the detuning is large enough. 

The above discussions on the discord dynamics inspire us to 
ask an interesting question: may the discord remain unchange 
in the time evolution before the critic time for certain cases? 
In the following we show that it is possible to maintain the 
discord before the critic time. In order to see this, we consider 
two different qubits with a large detuning r » 1, which are 
initially prepared X-type states with the three state parameters 
c\ - \, and < -C2 = C3 = c < 1. In this case, at zero 
temperature we find the critic time is given by 



f = _ Vc"2/(^^i)_i 



U)c 



(50) 



and the quantum discord has the following expressions 



T)(p) 



lf\og^{l+c)+'-f\og^{\-c\ t<tc, 

i^ iog2(i + n) + %i iog2(i - ri), f > fo 



(51) 



where y\{f) is given by Eq. ( |43] ). From Eq. ( fSTT l we can see 
that before the critic time the quantum discord is indepen- 
dent of the time. This indicates that the initial discord of the 
X-type states under present consideration does not change in 
the time evolution before the critic time tc- Hence, the en- 
vironment does not affect the quantum discord in this time 
evolution regime. In other words, the quantum discord is 
decoherence-free for such X-type states with the three state 
parameters c\ - I, and < -C2 = C3 = c < 1 in the time evo- 
lution before the critic time. From Eq. ( fSOl l we can see that 
the smaller the value of the state parameter c the longer the 
discord can remain unchanged. 



V. CONCLUDING REMARKS 

In conclusion, we have studied the quantum discord dynam- 
ics of two uncoupled qubits immersed in a common heat bath. 
This system is depicted by an exactly solvable phase deco- 
herence model. We have shown that the quantum discord of 
two noninteracting qubits can be amplified or protected for 
certain initially prepared X-type states in the time evolution. 
Especially, it has been found that there does exist the stable 
amplification of the quantum discord for certain X-type states 



8 



when the two qubits are identical, and the quantum discord 
can be protected for the case of two different qubits with a 
large detuning. However, the stable amplification cannot be 
created when the two qubits are not identical for the same X- 
type states. The degree of the discord amplification decreases 
with the increase of the detuning between the two qubits. The 
amplification and protection of the quantum discord can be 
understood as quantum effects induced by the environment 
since the two qubits interact neither directly , nor with a third 
system, but with the common reservoir with infinitely many 
degrees of freedoms. Both of them are decoherence free. The 
stable amplification (protection) of the quantum discord can 
also be considered as a resonant (large-detuning) effect of the 
two qubits interacting with the common reservoir 

It has been indicated that in general there does exist a sud- 
den change of the quantum discord for the two qubits in the 
time evolution at a critic time point tc, and the discord am- 
plification and protection may occur only in the time interval 



< t < tc for certain X-type states. Generally, the critic time 
depends on the initial-state parameters, characteristic param- 
eters of the two qubits and their environment, and the qubit- 
reservoir interaction. The quantum discord will decay asymp- 
totically to zero after the critic time. The creation of the stable 
amplification and protection of the quantum discord for the 
two qubits in a heat bath sheds new light on production of 
quantum states with long-living quantum discord. 
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